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Abstract 

We have proposed, in our previous papers]^, |2[, a method to characterize inte- 
grable discrete sohton equations. In this paper we generahze the method further and 
obtain a g-difference Toda equation, from which we can derive various g-difference 
sohton equations by reductions. 

1 Introduction 

It has been known that, for a given soliton equation, there exists a discrete analogue of 
the evolution equation, which preserves integrability. Once we find a discrete integrable 
equation, we can derive infinitely many integrable differential equations by taking contin- 
uous limits of variables in many different ways. We are interested in characterizing such 
integrable discrete systems. 

In our previous papers || we studied various integrable discrete systems which have 
certain symmetry and satisfy periodic boundary conditions. In this approach the nature 
of a quadratic equation, imposed by the boundary condition, plays a crucial role. The map 
which generates a time evolution turns out to be twofold, corresponding to two solutions 
of the quadratic equation. We showed that the discriminant of the quadratic equation 
becomes a perfect square when the system is integrable. Owing to this fact the map is 
free from a square root, hence is rational. We will call this type of map a non square root 
map (NSRM) in what follows. The discrete versions of Lotka-Volterra equation, KdV 
equation, Toda lattice, KP equation and Painleve equations 0, ^, ^, ||, ^ |, ||, |10[ are 
such examples. The purpose of this paper is to generalize this scheme of characterizing 
discrete integrable systems. 

To clarify the point of our argument let us show first a typical example. The 3-point 
discrete time Toda lattice is given by the set of coupled equations: 



XkUk — UkXk-l 



k = 1,2,3. (1) 



Here x^ means the variable Xk at the time t -|- 1, i.e., Xk{t) = Xk{t + 1). We impose the 
periodic conditions Xk+3 = Xk and Uk+3 = Uk. Hence there are 6 dependent variables. 
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Before solving eq.(|I]) directly it is more instructive to solve 



XkUk = Pk j^2) 

Xk + Uk+l = Qk 



for arbitrary functions pk, Qk of {x,u). We can also write them as 

QkXk+l — Pk+l 



Xk 



Xk+i 



Qk-iUk-i - Pk-l 
Uk = z • 

Uk-l 

Note that they are special cases of Mobius transformation. Hence by repeating this 
transformation, say for the first equation of eq.(^), three times, Xk is given in terms of 
(xfc+3, Uk-s), which are nothing but {xk, Uk) themselves owing to the periodic boundary 
conditions. In this way we obtain an equation 

Axk + B 



or, equivalent ly, 
where 



rxi -{A- A)xk -5 = 0, (5) 



A = qiq2qz - QkPk - qk~iPk-i, 

B = -QkQk+iPk+i + Pk+iPk-i, 

B = qk+iQk-i — Pk', 

A = -Qk+lPk+l, 



in the case of three point Toda lattice. 

The quadratic equation eq.@ has two solutions 



= ^ , (6) 

unless its discriminant 

Bis = {A- Af + ABB (7) 

vanishes. The general form of solutions contain square root of polynomial functions. 
Under this circumstance a sequence of the map will yield very complicated orbits in 
general. But the integrable systems are not such cases, i.e., the discriminant happens to 
be a perfect square of some polynomial. In fact if we substitute the right hand sides of 
eq.(|l]) into p^'s and g^'s of the expression eq.(|^), the discriminant is considerably simplified 
and is given by 

Bis = (X1X2X3 — U1U2U3Y. (8) 

This is a typical non square root map (NSRM), which we are going to study (see below 
for more precise definition). 
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This remarkable feature could be understood from other point of view as well. Namely 
by an inspection of our original equations eq. (|I|) we can see easily that 

Xk = Uk, Uk = Xk-1 (9) 

is one of sets of solution. This implies that another set of solutions of the quadratic 
equation eq.(^ must not have a square root either. We would like to emphasize that the 
existence of the solution of such simple form eq. @ owes to the symmetry of the equation 
eq.®. 

We will study, in the following section, A^-point maps under which the square roots 
are removed when periodic boundary conditions are imposed. In section 3 the discrete 
time Toda equation is studied in detail based on the preceding argument. In particular we 
will show that this scheme enables one to generalize the Toda equation to a g-difference 
form, so that various g-difference soliton equations are derived by reductions, including 
the qPiv equation by Kajiwara, Noumi and Yamada0. In section 4, we will discuss time 
evolution of solutions. The last section is devoted to summarize arguments in this paper. 



2 A^-point Map with Periodic Boundary Conditions 

We consider A^-point systems in this section and study conditions under which a map 
becomes a non square root map (NSRM). 



2.1 Mobius transformation 

Generalizing eq.@, let us consider the following sequence of Mobius transformations, 
Xk+hin-i) = "I" n = l,2,---,N, h = ±l, anSn-f3nln^0 (10) 

JnXk+hn + On 

and Mobius matrices 



M„ = ( ^" ) , det M„ 7^ 0, n = 1, 2, ■ ■ ■ , iV. (11) 



7n 5; , 

After times of the transformation we obtain the equation 



_ _ Axk+hN + B ^^2^ 

FXk+hN + ^ 

while the corresponding Mobius map is given by 

5^:= ^ ^ ) =MiM2---Mm. (13) 

We now impose the periodic boundary conditions Xk = Xk+hN- Then the formula 
eq.([T2D is equal to the quadratic equation 

rxl -{A- A)xk -B = 0. (14) 
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The solutions of this equation are given by 

ti Sn -2A±^/Drs 
= ^ (15) 

whereas the discriminant Dis is 

Dts = {A- Z\)2 + ABr = (tr Sn^ - 4 det Sn- (16) 

2.2 Non square root maps 

In this subsection we study maps such that the discriminant eq.(0) becomes a perfect 
square of some polynomial function. As we will prove shortly the following 6 patterns of 
the map satisfy the requirement: 

( Pn\ f an \ ( an Pn\ 

[ (3n o r [ 6n r I 6n r 



an \ -fn + (3n -/3n\ /3„ 

7n 5n / ' V 7n / ' I -7^ /?„ + 7n 



(17) 



Proof: 

Since the trace and the determinant are invariant under the transposition ot the ma- 
trices, we do not consider here the fourth and sixth patterns. For the first three patterns 
we find 

where 

a = ai---a]y, j3 = j3i- ■ ■ 5 = 5i---5n-, 
and F is a certain function of 5„. From this expression we easily find 

Dis = {2(3f5NMd. (a - ^f. (" - 

corresponding to the first, second and the third patterns. 
For the 5th pattern we first observe 

det Sn = det Mi ■ ■ ■ det Mn = (3^. (19) 

In order to calculate trS'Ar, we split M„ into two parts 



^'^"=( o" j + (;: oj-«»+^- (20) 

Using the property 

BmCn = 0, m,nel,2,---N, 
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we can see that Sn has terms of the form Ci - ■ ■ CpBpj^i ■ ■ ■ Bn only. From this fact, the 
trace of Sn is calculated as 

tr^^ = 5^tr(Ci---C,5,+i-.-5^) (21) 
= ^tr(C2---C,fi,+i---5^Ci) 

= tr(Bi---B^)+tr(Ci---C,v) 

= /3 + 7 



where 7 = 7i ■ ■ ■7Ar. Similarly we can treat the 6th pattern and obtain the same results 
as eqs.(|19|) and (^) for the determinant and trace, respectively. 

Hence we have obtained discriminants for all six patterns as follows: 



1) 


tr Sn 


^ N ,evenl3 1 


det Sn 


= -/5^ - 


L)zs = 


{2j3y6N,odd 


2) 


tr Sn 


= a + 5, 


det Sn 


= a5 


Dzs = 




3) 


tr Sn 


= a + 5, 


det Sn 


= a5 


Dis = 


(«-^)^ 
(7-/?)^ 


4) 


tr Sn 


= a + (5, 


det Sn 


= a5 


Dis = 


5) 


tr Sn 


= 7 + A 


det Sn 


= /?7 - 


Dis = 


6) 


tr Sn 


= 7 + A 


det Sn 


= /?7 - 


Dis = 



Therefore every discrete equation with the Mobius matrix pattern of eq.(|T^ has no 
square root, although they are not all. 



2.3 Rational maps 

Solutions of the quadratic equation eq.(|l4D are given by eq.(|l5|). The purpose of this 
subsection is to write them explicitly in the case of NSRM. First we note that the 2nd 
and the 3rd pattern as well as even case of the 1st pattern do not form a quadratic 
equation due to -T = 0. Hence it is enough to consider the odd case of 1st pattern, 4th, 
5th and 6th patterns. 



2.3.1 Odd N 1st pattern 

The odd case is given by the following setup: 

= ° Q ) , tr 5^ = 0, det 5^ = Dis = {2(3^. (22) 

Hence the solutions of the quadratic equation are 



_ ti Sn-2A±VdIs ±2P , , 

Xk = ^ = TT = =tl- 23 

2r 2/3 ^ ^ 
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2.3.2 4th pattern 

In this pattern we have 

S'at = ^ ^ ^ ) ' ^''^ = C( + det Sn = C(6, Dis = {a — 5)^. (24) 
To calculate F, we split M„ to two parts An and Dn as follows: 

Using the fact AmDn = 0, we see that Sn has terms of the form Di ■ ■ ■ Dp^iAp ■ ■ ■ Afq 
only, where p = 1,2, ■■■,iV + l and Dq = Aj^+i = 0. All possible terms can be written in 
such forms as 

A,A,---A^ = (^^ (26) 
DiD2---Dn = ( . ° ° ) , (27) 

Di---Dp^iAp---AM = [ ° ° I . (28) 

Hence we find the expression of F as 

N 

F = ^5i---(5p_27p-ittp---a7v- (29) 
The solutions of the quadratic equation are then given as follows: 

a — 5 

tr ^jv - 2Z\ ± \^Dis (a - 5) ± (a - 5) I , , 

= ^ = ^ = <i . (30) 







2.3.3 5th pattern 



We already got tr Sn and det Sn in the previous subsection, but did not have calculated 
the elements of Sn- To obtain the solution of the quadratic equation, we must find their 
elements since information about C and D are necessary for the solutions. 

The splitting idea of eq.(^) is convenient to obtain the elements here again. Using 
this idea , we need to calculate only the following quantities 

C,C2---Cn = l(] n 1, (31) 



1 

1 - 




B,B2---Bn = f^i I ^M, (32) 



Ci - ■ -Cp-iBp- ■ ■ En = 7i • ■ ■ 7p-i/9p ■ ■ ■ /?Af ( | _| ) ' (33) 
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because BmCn = 0. Hence F and A of the elements of Sn are obtained as 

N 

r = 7i---7iv + ^7i72---7p-i/5p---/5jv, (34) 

p=2 

N 

^ = -E7i72---7p-i/5p---/5iv. (35) 

p=2 

Consequently, noting the fact F — j — A, we obtain the solutions of the quadratic equation 
as 

tr -Sjv - 2Z\ ± _ (7 + /9) - 2/1 ± (7 - /3) 

- ^ - 2(7 - A) 

1 

P-A ^ PN QNiPn) ^^^^ 

[ -f-A - 71 Qi(/3,7) 
where Qn{P, 7) is the homogeneous polynomial of (/3i, • • • , /3jv, 71, • • • , 7Ar), that is 

N 

QniP, 7) = QniPl, • • • , /^AT, 71, • • • , 7iv) = X] 7n+l ' ' ' In+p-lPn+p ' ' ' (^n+N-1- (37) 

2.3.4 6th pattern 

For this pattern too, we already got tr Sjs and det in the previous subsection, but did 
not have calculated the elements of Sjq. Using the splitting idea again, the elements of 
S'at in this pattern are able to be obtained. The splitting of Mn in this case is 

Mn-i ' l = fn ^"l + f ' M-^n + C.. (38) 

\-ln Pn + ln J \ Pn J \ 'In In ) ^ ' 

Since C^Bn = 0, we are enough to calculate the following terms 

CiC2---C^ = 7(_°i 5), (39) 
B,B,..-B^ = /?(^° \y (40) 

B,---Bp^,C,---Cn = PiP2---Pp-i1p---1n( Z\ []■ (41) 



Hence we obtain F and A of the elements of -S'jv as 



N 

r = -7 - I] A • • • Pp-ilp ■■•In, 

p=2 

N 

A = i3 + ^+Y.P^---^-^^p---^N. (42) 

p=2 



If we further use F = P — A, we get solutions of the quadratic equation as follows: 

trSN -2A±VDi^ (7 + /3) - 2Z\ ± (/3 - 7) 



2r 



2{f3-A) 



1 

Z\ - 7 
A^ 



1nQn{i,(^)' 



(43) 



2.4 Polynomial Q 



We note that polynomial Q's, which appear in the solutions of the 5th and 6th patterns, 
are given by the minor determinants of the diagonal entries of the following matrix. 



<5(/5i,---,/3jv,7i,---,7iv) := 



/A +71 A 

72 (^2 +72 (^2 

73 /^S + 73 








7i 















JN Pn + J 



For simplicity, we will introduce the following notations. 



(44) 



Q^(/9,7) 



Q{pi,---,pN,ji,- 
Q(7i,---,7iv,/3i,- 

Q{f3N: ■ • •,/5i,7Ar, 

QhN,---m,pN,- 



■,1n), 

■,7i), 
•,/9i). 



Namely, the matrix Q having entries in ascending order is expressed by Q, and in de- 
scending order by Q^. These matrices are related by 



where 



g^(/3,7) = jg(7,/3)^, 
/ ••• 

J := : 

1 

V 1 



1 \ 

1 



(45) 



J 
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We set symbols of them of (n, n)-entry as Qn, ^-e. 

N p-1 N-l 

N p-1 N-l 

p=l p=l i^=p 
N p-1 N-l 

Qnif^n) = H n In+N-t, n Pn+N-u, 

p=l fj,= l v=p 

N p-1 N-l 

"Jn+N-ui 

p=l p=l y=p 

where j3n±N = Pn, 7n±N = In- 

These polynomials satisfy the following useful relations. 

1. QN-n+l{P^l)=Ql{l.(3)^ Qn(A7)=Q)^-„+l(7,/?), 

2. {Pn + ln)Qn{P,l) = f3n-lQn-l{f3 , l) + ln+lQn+l{f3 , l) ■ 

The relation |I] is obvious. The relation ^ is non-trivial and will be checked by means of 
more fundamental relation : 



PnQn{P,7) = /?-7 + 7n+l<5n+l(/3,7) 

where (3 = (3i - ■ ■ (3n, 7 = 7i ■ ■ ■ 7Af- 

We can show that Qn{P,'y) coincides with the cofactor of the (n, n) element of the 
product RL of the following matrices 



/ 1 
Pi 1 
P2 I ■■■ 



Pn \ 



/ 7i 1 • • • \ 

72 1 



R 



V ■■■ Pn-i 1 / 
which could represent the dToda equations 





V 1 



In J 



2.5 Modification 

The discriminant Dis is invariant under arbitrary similarity transformations, since it is 
given by the trace and determinant of S^. This fact enables us to modify our map starting 
from the six patterns we considered above. 
Namely under the transformation of M„, 

M„ ^ = U;,^Mr,Un+i, (46) 

the matrix eq.([TB|) is changed into 

^^"^ = Mf'^Mf"'^ ■ ■ ■ M'^"'^ = t/f ^MiMa ■ ■ ■ M^Un+i • (47) 
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If we take into account the periodic boundary condition, we must have Un+i = Ui. 

Therefore the trace and the determinant of Sn, hence the discriminant Dis as well, 
remain constant: 

tiS'^"'^ = tTSN, detS^""^ = detSN, Dis"""^ = Dis, (48) 

under the transformation eq.(^) for arbitrary matrices UnS. This implies that the mod- 
ified map M™°'^ is NSRM if M„ is so. 

Now denoting by x"^°'^ the modified variable of a;, we will write the Mobius map 
associated with M„ as Mn{xk)- In this way, the modified Mobius map is possible to be 
written as 

^w^-mod^ = ?7f ^Mit/i(a;^+°f). (49) 



X 



mod 
k 



Here if we define = U^li{xk+n)^ then 



T^mod 



X 



U{'M,U2{xff,) 



Ur\xk) = U^'MiU2U^\xu+i) (50) 
Xk = Mi(xfe+i). 



As a result, the following two modifications are equivalent 



2) Xk+n ^It+n ~ ^n+l(^fc+n) 



This correspondence is preserved in the solutions of the quadratic equation, so we can 
create many equations from the basic six patterns in this way. 



3 Discrete Time Toda Equation 

We study, in this section, the A^-point discrete time Toda equation (dToda) in detail from 
the view point of the previous section. 



3.1 Two types of solution 

Let us start with writing the A^-point dToda equation. 

I ^'^'''^ " ''^^''-^ , A; = 1,2,3, AT. (52) 

[ Xk + Uk+l = Uk + Xk 

In order to apply the argument in the previous section to this system we transform eq. 
into the form of Mobius transformation. 



Xk + Uk)Xk+l — XkUk+l 




Xk+1_ 

Xk-l + Uk-l)Uk-l - Xk-2Uk-l ^ ' 

Uk-l 

Xk,Uk). 
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We notice that this is a modified 5th pattern Mobius transformation which is specified 
by the following data: 



+n— Ij Uk-\-n—l ) 
{Uk—ni Xk—n) 



1 

7n 



(54) 



Substituting these data into eq.(^) we obtain solutions immediately. Before writing 
them down, however, it is convenient to simplify notations. Namely the matrix Q, which 
appears in the solutions, is Q{xk, Xk+i, ■ ■ ■ , Xk-2, Xk-i, Uk, Wfc+i, ■ ■ ■ , Mfc-2, Wfc-i)- Let the 
matrix Q be 

Q{Xk,Uk) := Q{Xk, ■ ' ' :Xk+N-l,Uk, ■ ■ ■ ,Uk+N~-l), 
Q^{Xk-,Uk) := Q{Xk+N~l-,- ■ ■ iXkiUk+N-l-,- ' ' i^k)-, 

and its minor determinants be 



Qn{Xk,Uk) :— QniXk, ■ ■ ■ ,Xk+N-l,Uk, ■ ■ ■ ,Uk+N-l), 
QniXk,Uk) := Qn{Xk+N~l, ■ ' ' ^ Xk,Uk+N-ir ' ' jUk). 

These minor determinants satisfy 

1- QN-n+l{Xk,Uk) = QniUk, Xk), QniXk, Uk) = Q'if_n+i{Uk, Xk) , 

2. {Xk+n + Uk+n)Qn{Xk, ttfc) = Xk+n-lQn-l{Xk, U^) + Uk+n+lQn+l{Xk, Uk), 

3. Qn{Xk+N,Uk+N) = Qn+l{Xk+N-l,Uk+N-l) = ■■■ = Qn+N-l{Xk+l,Uk+l] 
Qn+N^Xk, Uk) 

where Xk±N = Xk, Uk±N = Uk- Therefore the solutions turn out to be 

QN{Xk,Uk) 

A-type, 



(55) 



Xk Xk—1 



Uk = Uk 



QliXk,Uk) 

Ql{Xk,Uk) 
QN{Xk,Uk) 



(56) 



Xk = Uk 
Uk = Xk-1 



B-type. 



Similarly the time reversal map is obtained by writing eq.(0) as 

{Xk + Uk+l)Xk-l - XkUk 



Xk 
Uk 



Xk-1 

{Xk + Uk+l)Uk+l - Xk+lUk+1 
Uk+1 



(57) 



corresponding to the map {xk-i,Uk+i) — > {xk,Uk)- 

This is again a modified 5th pattern of the Mobius transformation. The modification 
matrix and the elements of the map are given by 



{Xk+l~n, Uk+2-n) 
{Uk+n, Xk-l+n) 



1 

7n 



(5^ 
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and the solutions 



•^fc— 1 Xk 



Ql{Xk,Uk+l) 

QX{xk,Uk+i) 
Q^{xk,Uk+i) 



Xk — Uk+l 
Uk = Xk 



A-type, 



B-type. 



(59) 



Instead of eq.(|53|) and eq.(|5^) we could write the dToda equation eq.(|5^) in the form 

Xk -- 



Xk 
Uk 



-Xk-lUk 


Xk-1 


- (Xfc-l + Uk-l) 




Xk-lUk 


Uk+1 


- {Xk + Uk) 




—Xk+lUk+1 


Xk+1 


- {Xk+l + Uk+l) 




— XkUk 




- {Xk-1 + Uk) 


and 


{Xk+l,Uk-l) 



(60) 



(61) 



6th patterns whose elements and the modification matrix are given by 



(/3n,7n) 



(•^fc— 1—715 U,k—l~n) 
{Uk+ nj Xk+nj 
(^^fc+l+ji; Uk+l+n) 
{Uk—n^ Xk~l—n) 



Un 



l/ln-l 

1 



(62) 



for eq.(pD|) and eq.(^), respectively. Applying these data to eq.(^3|) we obtain the same 
results of eq. (^) and eq. (^9]) . 



3.2 Generalizations 

3.2.1 Symmetry of B-type solutions 

We have seen that the B-type solution is simply 1 if the map is the form of either 
5th or 6th pattern. A modification of the matrix by f/„ will change the solutions as well 
as the equation while the discriminant remains constant. In particular a different B-type 
solution will correspond to a different set of equations, and hence different A-type solution. 
If the B-type solution is rational the A-type solution must be also rational satisfying the 
same equations. 

To see how it works we try to generalize the B-type solution eq. (^6]) to 



Xk — Uk+i 



I Uk Xk+j 

It is not difficult to find their corresponding equations 



0,l,2,---,Ar-l. 



Xk+lUk — Xk+jUk+i+l 
Xk ~l~ Uik^m Xk+j+m 



Uk+i 



I, m 



0,l,2,---,Ar-l. 



(63) 



(64) 



12 



From the construction this equation has also A-type solution which is a rational polyno- 
mial of X and u. And we can set h = m + I = ±1 without loss of generality. 

We would like to emphasize that the B-type solution follows automatically to the 
symmetry possessed by the equation. In other words, the rational solution of A-type is a 
consequence of the symmetry of the Toda lattice equation. 



3.2.2 g-difference version 

In the paper Kajiwara et al. have studied affine-Weyl symmetry of the discrete time 
Painleve IV equation with parameters, which is called qPjv- We have shown in our 
previous paper that gP/v is one of the equations which are characterized by the square 
root free maps. 

Now we want to know whether exist similar generalization of the discrete time Toda 
lattice eq. (^) which we have been studying. To this end we consider the following equa- 
tion with a set of parameters ctfc, 6^, = 1, 2, ■ ■ ■ , iV : 

T~ ~ Clk+jXk+jOk+i+lUk+i+l 

_ cik+iyk /•g5^ 

•^k . Uk-\-m \ h 

r — — flfc+j+m^^fc+j+m + ^k+iUk+i- 

O'k C'fc+m 

Let us call this equation a g-difference generalized discrete time Toda equation (gToda). 
Remarkably this equation has the symmetry of modified 5th or 6th pattern of the Mobius 
map. So, this is transformed to a quadratic equation of every x, u, x or u, and has two 
types of map as solutions to eq.(|65|). 
In fact if we introduce 

j (^k+s+nh-^k+s+nh ~l~ ^fc+j+n/i^fc+j+nh (^k+s+nh-^k+s+nh j (^66) 

Y bk+i+nh^k+i+nh J 

Un = (I I ], (67) 

\ U ^k+i+nhU-k+i+nh ) 

we obtain two maps as follows : 



a-k+s-hXk+s-h-p^, 7 T, A-type 

Vl [O'k+sXk+si Ok+iUk^i) 



(6^ 



^ bk+iUk+i, B-type 
where s = j + m,h = l + m, and is defined by 



Qn{xk,Uk) :-- 



Qn{Xk,Uk) if h=l 
Qn{Xk,Uk) if h = - 



3.2.3 Higher dimension 

We can further generalize gToda to higher dimensions by increasing the number of suffices: 

{Xk,Uk) (yXk-i,k2,---,kMi Uki,k2,---,kM) ■ 
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The B-type map should have the form 



(^ki,k2,---,kM 
'^ki,k2,---,kM 



Jki,k2,---,kM 

In the case of M = 2, for instance, 

r Xk£ 



Uk,l 

hhi 



bki+ii,k2+i2,---,kM+iM'^k-i+ii,k2+i2,---,kM+iM 



(^ki+ji,k2+j2,---,kM+jM^ki+ji,k2+j2,---,kM+jM- 



^k+l,l+lUk+l,l+l 



O-k+ld-l^k+ld-l 



satisfies 



Xk-l,lUk-l,l 
dk-ldbk-lA 



Xk,l _|_ Uk-l,l+l 



0-k,l-lbk,l+lXk,l-lUk,l+l, 



. , 0'k,lXk,l + bk+l,l+lUk+l,l+l- 

. (lk,l Ofc-l,/+l 

This is the gKP equation which we derived in our previous paper. 



(69) 



(70) 



(71) 



3.3 Reductions 

It has been well known that the discrete time Lotka-Volterra equation (dLV) can be 
obtained from the discrete time Toda equation by a Miura transformation |^ . We like to 
show, in this section, that the same dLV-type as well as the discrete time KdV (dKdV)- 
type equations can be derived from the discrete time Toda lattice eq.(^) by simple 
reductions. The basic idea of our reduction is to deal with one of two equations in gToda 
constraint to eliminate half of the variables. 



3.3.1 Reduction to the dLV-type 

The way of reduction to a dLV-type equation is to impose an additional constraint to the 
second equation of gToda as follows. 



Xk+l Uk 
O-k+l bk 



(^k+j^k+jbk+i+lUk+i+l 

(72) 



■^k Uk+m 7 

r 7 = ak+j+mXk+j+m + Ok+iUk+i = — Tfc 

O'k bk+m 

where r^'s are constants. 

This constraint enables us to eliminate variables, say u and u, from the first equation. 
Namely the following substitution 

Uk+m ^ -^k ^ 

bk+m O-k (73) 

bk+i^k+i ^ (J'k+j+m-^k+j+m ^k 
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yields the dLV-type equation 



■^k+l I •^k—m 
I ^k—m ~r 



O-k+l 



a-k- 



O'k+jXk+j [f^k+l + (^k+j+m+lXk+j+m+l) ■ 



(74) 



When rk = ttk = bk = I, k = 1,2, ■ ■ ■ , N, this is the dLV equation. Various types of 
gLV equations which were discussed in our previous paper, are included in this equation. 
It is straightforward to obtain two maps of this equation. We find 



Xk_ 

dk 



Q%{(^k+s^k+si ^k (^k+j+m^k+j 

^k+s—h-^k+s—h 



A-type 



Q\i^Ojk+s-^k+si Ojk+j+m-^k+j+va) 
1"k ~ 0,k+j+m.Xk+j+m-, B-type 



(75) 



which could be also obtained from the map of qToda eq. (|65|) by the substitution eq. (|73D . 
The gPiv of Kajiwara, Noumi, Yamada is included as one of the A-type map. 

3.3.2 Reduction to the dKdV-type 

Let us try an alternative case of constraint. Namely we impose a constraint to the first 
equation of eq.( 



Xk+l Uk 
dk+l bk 
■^k ^ '^k+m 
O'k bk+m 



O'k+jXk+jbk+i+lUk+i+l 
O'k+j+m-^k+j+m ~l~ bk-i-iUk+i- 



rk, 



(76) 



This enables us to eliminate u and u from the second equation according to the rule 

O'k+l'f' k 



Uk 

bk 

bk+i+lUk+i+l 



Xk+l 

rk 



(77) 



(^k+j-^k+j 



from which we obtain 



•^k , (^k+l+mrk+m 

\ = = flfc- 

O'k Xk+l+m 



-j+rrt 



■^k+j+m ~l~ 



Tk-l 



dk+j-lXk+j- 



(7J 



This equation containes the discrete time KdV equation as a special case. Hence we 
call eq. (|78D the gKdV equation. The corresponding solutions are given again by simply 
replacing u and u in the solutions of gToda according to eq. (ff^) . 



Xk_ 

dk 



^k+s—h-^k+s—h ^fj^ 



Q^{,0,k+sXk+S) fk-l/dk+j-lXk+j-l) 



A-type 



(^k+s-^k+sj ^k— l/ O-k+j-lXk+j-l) 

, Tk-i/ak+j-iXk+j-i. B-type 



(79) 
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3.4 Reflection maps 



We define some symbols of the two types of map and their inverse maps of dToda as 
follows, 



A 




( QN{Xk,Uk) \ 



A- 




Uk 



Xk 



QliXk,Uk) 

QljXk, Uk) 
QN{Xk,Uk) J 

Q)^{xk,Uk+i) \ 
QX{xk,Uk+i) 
QX{xk,Uk+i) 



B 





B- 





' Q)q{xk,Uk+i) J 
These maps satisfy the following relations. 

(I) AA-^ = A-^A = 1 

(II) AB^ = B^A, 

(III) ABA = B'\ 



BB- 



B-^B = 1, 



or 



(I') 
(II') 

(III') 



BB-^ = B-^B 
CB = BV, 
C^ = V^ = 1, 



^0) 



(81) 



where C = A^^B = B^^A and V = AB^^ = BA^^. The relations eq. (^ll) are derived 
from eq.(^) by some calculations, automatically. 

The reflection maps C, T> change the inverse A-type map to the inverse B-type one 
and the A-type one to the B-type one, respectively. Namely, 



C : inverse A-type 
inverse B-type 

V : A-type 
B-type 



inverse B-type, 
inverse A-type, 
B-type, 
A-type. 



Hence, these reflection maps are permutations of solutions of the quadratic equation 
eq.(|l^ in the case of dToda, i.e. Galois tranformations of eq. (p!^ . Furthermore, the 
generalized version of dToda eq. (|6^) satisfy these relations eq. (|80D or eq . (^TD , but gToda 
does not. These reflection maps play an important role in affine-Wyel symmetries and 
g-difference Painleve equations as studied by Kajiwara, et al. p^ . 



4 Discrete Time Evolution of the gToda 

A sequence of B-type map does not generate an orbit but causes exchange of the variables x 
and u, whereas an A-type map generates an orbit. We will study a discrete time evolution 
of A-type map in this section. To make the time dependence explicit we write the variables 
X and u as x*"*"^ and m*"*"^ in the following. 
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4.1 Time evolution of the A- type maps 

A repeated use of eq.(|5BD enables one to express x^"^^ and u^,^^ in terms of x°'s and m°'s 



as 



X 



t+i 



X 



k-l' 



Qi 



X 



k-t-l 



n 



L Qt-S ' 
Vfc-s 



U 



t+1 



U 



t Qk _ _ ,,0 TT ^ 



(82) 



Here we denote by the following polynomial of a; s and u s, 



N 



Qk ■= 



* 7/* ■ 



'^k+j-l^k+j 



■ X 



fc+Af-1- 



^3) 



If we substitute eq. (|8^) into the right hand side of eq. (|83D we will find a considerable 
cancellation of factors in the numerator and the denominator. The cancellation might be 
associated with a reduction of algebraic entropy |Tl|, 0. As a consequence we obtain a 
recurrence equation for Q*: 



N 



■ "^k+j-lXk+j-t-l 



' X I 



■j-j- Vfc+j-iVfc+j-s-i 

s=0 Qk Qk+N~s~l 



The A-type solutions are given if is found by solving eq.(|^) iteratively. The 
generalization of this result to the case of gToda is straightforward. This is, however, 
a cumbersome task. And so, we will look orbits of A-type map described by numerical 
computation. 



4.2 Numerical observation of ultra-discrete version of the A- 
type maps 



It is not difficult to see behaviour of orbits of the A-type map eq.(p2D as long as we generate 
them numerically by using computers. They are constrained on curves since there are 
sufficient number of conserved quantities. If parameters a^'s and bkS are introduced the 
orbits tend to expand and closed curves turn to spirals. These are behaviours which we 
could expect naively from our previous experiences. Instead of studying them in detail 
we would like to present here behaviour of ultra-discrete version of solutions. 

It has been known that if a discrete integrable equation has its ultra-discrete version 
the integrability is preserved |T^. There exists a systematic method to ultra-discretize a 



discrete integrable equation if it contains only product and summation of variables but 
not subtractions. Namely we introduce a new variable X and replace x by a; = e^^^. The 



ultra-discretization procedure [T^ is defined by 



lim (e^/^ ■ e^/') ^ X + Y, lim (e^/^ + e^^') max(X, Y). 

To be specific let us consider the case of dKdV-type equation. The g-difference version 
has been given by eq. ([78|) . If we further specify to 3-point case the equations are 

Xk ^fc+l 1 

h = cLk+i+jXk+i+j H , (85) 

O'k Xk+i ak+jXk-\-j 
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Figure 1: Ai = A2 = A3 = Figure 2: j = 0, Ai + A2 + A3 = 




fc — 1,2,3, j — 0,1,2, — Ofc, Xk+3 — Xk- 

The corresponding A-type map can be also found readily from the one of g-Toda equation 
by the reduction eq.(^). 

1 2 2 

1 + 0'k+20'k+lXk+2Xk+l + CtkO'k+l^k+2XkXk+l^k+2 

Xk-j — 0'k~jO'k+2Xk+2 — ; 2 2 ' v""/ 

i + akak+2XkXk+2 + ClkO'k+lO'k+2^kXk+lXk-^2 

k = 1,2,3, j = 0,1,2, ak+3 = ttk, Xk+3 = Xk- 
The ultra-discrete version of eq. (^6]) turns to 



^k-j — ^k-j + Ak+2 + Xk+2 

+ max(0, Ak+2 + ^fe+i + Xk+2 + Xk+i, Ak + 2Ak+i + + Xk + 2X^+1 + Xk+2) 
- max(0, Ak + Ak+2 + Xk + Xk+2, Ak + Ak+i + 2Ak+2 + Xk + Xk+i + 2X^+2). 

(87) 

Here we set ak = e"^^/"^, k = 1,2, 3. Let X^'s and AkS be rational numbers and impose 
a condition Ai + A2 + A3 = so that the map returns to initial values after some steps. 
Under these circumstances we find maps such as presented in Figures |l] to 

Figure 1 shows the case of = 0, k = 1,2,3 for all j=0,l,2. Similarly Fig. 2 
corresponds to the case of + ^2 + A3 = and j = 0. The straight lines in these figures 



are known to correspond to integrable orbits [jT4[. Figures 3 and 4 are the cases of j = 1 
and j = 2 respectively. They do not seem to describe integrable orbits. Therefore the 
orbits behave differently for different value of j. This must be contrasted with the case 
of dKdV in which orbits form smooth closed curves under the condition 010203 = 1 for 
all j = 0,1, 2. 

It will be worth while to point out that the dKdV map eq. (^6]) in the case of j = is 
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Figure 3: j = 1, Ai + A2 + A3 = 



Figure 4: j = 2, Ai + A2 + A3 = 




invariant under the following Backlund transformation 

7r(afc) = Ofe+i, TT{xk) = Xk+i, 



Si[aj) = aj[ai+iai^iy"\ Si[Xj) = x 



7ii 



1 + aj_|_iaj_iXj+iXj_i 



where 

2 -1 -1 \ / 1 

A := (A.,) =1-1 2 -1 , r := (7,,) = -1 

-1-12; V 1 -1 

and A; = 1,2,3 and ak+3 = ctfc, Xk+3 = Xk- This is the same affine Weyl symmetry- 
shared by qPiv 0- But A^(>4)-point dKdV maps of j = does not have this symmetry 
even though it exhibits an integrable orbit. 

We have also studied 3-point ultra-discrete LV-type maps in the same manner. The 
results are similar to those of KdV-type maps and found one integrable map corresponding 
to the qPiv- 



5 Summary 

The main purpose of this paper was to clarify the method characterizing discrete integrable 
systems, which we introduced in our porevious works. We have shown that, in the case 
of integrable systems, two adjacent variables are related by a Mobius transformation of 
special forms. We found at least 6 patterns of such forms. 

This observation enabled us to obtain quite large class of new discrete integrable sys- 
tems. The generalized discrete Toda equations (dToda) of eq. (^) is one of such examples. 
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By certain reductions of the dToda equation dKdV and dLV were obtained. The method 
admits introduction of arbitrary parameters into the equations and further generahzation 
to gToda systems. 

If we solve the equations in the form in which each variable is given as a function of 
variables in the previous time, there exist always two types of solutions which we called 
A-type and B-type. The B-type map does not create orbit but exchanges variables, while 
a sequence of the A-type map generates an orbit. 

The map of A-type is written as a rational polynomial of variables whose numerator 
and denominator are expressed by minor determinants of the matrix Q in eg . (|4^) . It is 
interesting to notice that this matrix Q is the same one known in the field of affine Weyl 
group and g-Painleve equations 0, |10[ . 

We have not succeeded analyzing full behaviour of a sequence of A-type maps. It 
is one of the most interesting problems to find behaviour of the map when parameters 
Ofc's and fo^'s are introduced. Some numerical results using ultra-discrete method have 
shown that the introduction of the parameters preserve integrability in some cases but 
not always. This result seems to show that the non square root map condition is not 
sufficient but necessary for a discrete map being integrable. It is desirable to clarify this 
point further. 
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